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ABSTRACT 

The structure of unstable planetary waves is computed by a quasi-geo- 
strophic model extending from the surface up to 80 km by means of eigenvalue- 
eigenfunction techniques in spherical coordinates. Three kinds of unstable 
modes of distinct phase speeds and vertical structures are identified in the 
winter climate state: (1) the deep Green mode with its maximum amplitude in 

the stratosphere, (2) the deep Charney mode with its maximum amplitude in the 
troposphere', and (3) thq shallow Charney mode which is largely confined to the 
troposphere. Both the Green mode and the deep Charney mode are characterized 
by very slow phase speeds. They are mainly supported by upward wave energy 
fluxes, but the local baroclinic energy conversion within the stratosphere also 
contributes in supporting these deep modes. The mesosphere and the troposphere 
are dynamically independent in the summer season decoupled by the deep strato- 
spheric easterly. The summer mesosphere supports the easterly unstable waves 
1-4. Waves 3 and 4 are identified with the observed mesospheric 2-day wave and 
1.7-day wave, respectively. 

INTRODUCTION 

Since the discovery of the spectacular natural phenomena of stratospheric 
sudden warmings, interest in planetary waves has been increasingly enhanced. 
After satellite observations became possible, along with previous observations 
by meteor radar, partial-reflection radar and rocket sounding, many transient 
wave activities in the upper atmosphere have been documented, such as 2-day 
waves in the summer mesosphere (RODGERS and PRAIA, 1981), 4-day waves in the 
upper stratosphere (VENNE and STANFORD, 1982), the slowly moving planetary waves 
in the winter stratosphere and mesosphere (HARTMANN, 1976), and a wide wave 
spectrum in the lower stratosphere (YU et al . , 1983). Some wave activities, 
mainly in wave 1 and 2, are intimately related to the sudden warming phenomena, 
which extend from the stratosphere well into the mesosphere. The planetary 
scale waves possess deep vertical structure in winter, but vertically decoupled 
in summer (LAEITZKE, 1981a, b). 

The purpose of the present study is to investigate those transient plane- 
tary waves based on instability theory. The growth rate, phase speed, merid- 
ional and vertical structure of unstable modes are examined for the winter and 
summer solstice climate states. A comparison study of stabilities for the 
basic states before and after the 1976/1977 sudden warming event is also dis- 
cussed. 

METHOD OF ANALYSIS 


(a) Governing Equation 


Assuming the large-scale wave motions are adiabatic and nondissipative, we 
use the linearized quasi-geostrophic, spherical model formulated by MATSUNO 
(1971): 
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Equations 1 and 2 are vorticity and thermodynamic equations, respectively, where 
X denotes the longitude, 6 the latitude, c = cos 9, s = sin 0, ft the angular 
velocity of the earth, a the radius of the earth, 4> the perturbation geopo- 
tential, U the basic zonal flow, and S the static stability. 

(b) Grid Point Arrangement and Boundary Conditions 

The model is discretized vertically by ten computation levels covering 
from the surface up to 80 km and it is referred to as an MST model (M for meso- 
sphere, S for stratosphere and T for troposphere). The computation levels are 
arranged at p “ 0.01, 0.04, 0.2, 1 mb in the mesosphere; at p = 3, 10, 30, 100 
mb in the stratosphere and at p = 250, 750 mb in the troposphere. The vertical 
boundary condition to = 0 is assumed at p = 0 and p = 1000 mb. In the meridion- 
al direction the model is bounded by the fixed boundary conditions <j>' = u' * v' 

• 0 at the pole and the equator with 10° meridional mesh sizes. The grid point 
arrangement is shown in Figure 1 , 
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Figure 1. Vertical and meridional grid arrangement 
and boundary conditions used in the spherical model. 
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(c) Basic State Parameters 

The winter and summer solstice climate basic flows are shown in Figure 2 
The zonal-mean temperature field and the wind field consist of the parameter 
space (S, 0), where the static stability is calculated from the temperature 
field. 

(d) Normal Node Solution 
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We assume a normal mode solution: 
$ = 4(6. p )e im(X - Ct) 


and calculate the phase speed C , growth rate a. = mC., and the eigen- 
function $(0,p), where C = C + iC. and $ = $ + i$ . are complex 

numbers. * 

RESULTS 

(a) Winter 

Figure 3 shows the spectrum of the growth rates and phase speeds of the 
unstable waves computed for the winter basic flow. The values of a and a g 
are representative for "dry" and "wet" tropospheric conditions. Both reveal 
the presence of three kinds of unstable modes: 

1. Green mode, denoted by MST in Figure 3, for waves 1 and 2, 

2. Deep Chamey mode, denoted by ST in Figure 3, for wave 3, and 

3. Shallow Chamey mode, denoted by T in Figure 3, for waves 5-9. 



Figure 3. The growth rate (ordinate) and the phase speed at 
the equator (labeled in m/s) for winter solstice basic 
state: a) the dry troposphere; b) the moist troposphere. 
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The phase speeds given in Figure 3 reveal an interesting instability ex- 
change between the slow modes (waves 1-4, including the Green mode and the deep 
Charney modes) and the fast modes (waves 5-9, all of which are confined to the 
troposphere). It implies some intrinsic relation between the phase speed and 
the vertical structure. The vertical wave energy fluxes (Figure 4) depict that 
the deep waves are vertically coupled as an entirety. They may be identified 
with the slowly moving waves 1, 2 and 3 in the real atmosphere. 

(b) Summer 

The spectra of the growth rates, phase speeds for the summer circulation 
are given in Figure 5. The mesospheric modes, with large negative phase speeds 
C r < 0, and the tropospheric modes, with the positive phase speeds C r > 0, 
exchange their stabilities at wave number 4. The computed most unstable mode 
is wave 3 with a period of 2 days which compares favorably with the observed 
2-day wave. Figure 6 shows structures of the computed mesospheric waves with 
m “ 1-4 and the tropospheric wave with m = 6. They are spatially separated by 
the easterly wind in the stratosphere. The planetary scale waves dominate in 
the mesosphere and the synoptic scale waves dominate in the troposphere. 

VERTICAL COUPLING AND DECOUPLING 

In order to understand the stability due to vertical coupling between the 
mesosphere, the stratosphere and the troposphere in the winter season, a case 
study is performed to compare the eigenmodes in the basic states before and 
after the 1976/1977 sudden warming event. Figure 7 shows the observed varia- 
tion of geopotential amplitude for waves 1 and 2 during the warming phase. 
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Figure 4. Geopotential eigenfunctions (upper and vertical wave energy 
flux (lower) in the NH winter basic flow with the "dry" troposphere. 
Units are arbitrary. 
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Figure 5. The growth rates (ordinate), and the phase speeds 
(labeled in m/s) for summer solstice basic state with 
a) the "dry" troposphere and b) the 'Wist" troposphere, 
m indicates the wave number, M the mesospheric mode and 
T the tropospheric mode. 


Both waves 1 and 2 were growing during December 1976 and decreasing in January 
1977. The basic flow in December 1976 is characterized by a broad, weak 
throughout the troposphere and the stratosphere, and in January 1977 is 
characterized by easterlies at high latitudes. Figure 8 summarizes the results 
of frequency computations for December and January. 

Three significant changes occur after the sudden warming: 

1. The maximum growth rate shifts from planetary- scale (m = 2) to synop- 
tic scale (m = 7) and the planetary scale waves are substantially 
stabilized. 

2. Transition occurs from the slow modes to the fast modes at wave number 
4 before the warming, but all waves change to the fast modes after it. 

3. All the deep modes (MST modes and ST modes) are suppressed to the 
shallow, tropospheric modes after the sudden warming as shown in 
Figures 9 and 10. 

Since the deep Green mode and the deep Chamey mode have been calculated 
in the winter climate basic flow, we suspect that the deep modes exist in the 
whole season except in periods after sudden warmings, in which all deep modes 
are suppressed primarily by the reversed zonal flow at high latitudes. 
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Figure 9. Geopotential (top) and wav.e energy flux (bottom) 
in December 1976. Units are arbitrary. 









